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ABSTRACT 

We consider the polarization of 21cm line radiation from the epoch of reionization due to both 
intrinsically polarized emission and secondary mechanisms. We argue that Thomson scattering of the 
21cm quadrupole by the reionized universe is likely to produce the largest effect. The 21cm quadrupole 
is sourced by baryonic density fluctuations and the fluctuations in the ionization fraction due to dis- 
crete H II regions. Since Thomson scattering produces only E-type polarization and is achromatic, 
its unique fingerprint can in principle be separated from foregrounds associated with polarized syn- 
chrotron emission which should not be correlated with the cosmic signal. We estimate that Poisson 
fluctuations of H II regions at the end of reionization (zr ~ 6-20) produce a brightness temperature- 
polarization cross-correlation signal of ~ 0.1 — 0.3 mK on angular scales of tens of arcminutes. This 
cross-correlation signal is within the instrument sensitivities of the future Square Kilometer Array 
(SKA) and close to the sensitivities of the forthcoming Mileura Widefield Array (MWA) and Low 
Frequency Array (LOFAR). 

Subject headings: cosmology: theory - intergalactic medium - diffuse radiation 



1. INTRODUCTION 



Measurements of the power-spectrum of 21cm 
brightness fluctuations can potentially constrain 
the complex process es that r eionized the un iverse 
JSuuvaev fc ZeldoTich| 17973 [Horan fc Reed I7T)73 



Scott fc Reesl 1 1 !)!)(* iTozzi et alJ |2000t iMadaueta! 



19971 iFurlanetto et alJ 12004 IBarkana fc Loe 
2QQ§1j D as well as the ba s ic cosmological parameter s 
(|Loeb fc Zaldarriagal 12004" IBarkana fc Loebl l2005ald) . 

By tuning low-frequency radio arrays to different fre- 
quencies, one could observe the intergalactic medium 
(IGM) at different redshift slices as it evolved through 
the epoch of reionization. Observations of the tempera- 
ture brightness power spectrum will allow us to infer the 
baryonic power spectrum and the peculiar velocity field 
at high redshifts, as well as the mean cosmic ionization 
fraction and fluctuations in the ionization fraction due 
to discrete H II regions. Constraints on these quantities 
will test models of reionization, in particular the spectra 
and evolution of the fi rst sources of ionizing radiation 
ijBarkana fc Loebl 12001). Several low-frequency arrays 
are currently under construction for the purpose of map- 
ping neutral hydrogen during the epoch of reionization 
(http://space.mit.edu/eor-workshop/); these include 
the Frimeval Structure Telescope (PAST) , the Mileura 
Widefield Array (MWA) 5 , and the Low Frequency Array 
(LOFAR) 6 , that would establish the grounds for the 
construction of the Square Kilometer Array (SKA) 7 . 

Almost all theoretical work thus far has focused on 
the brightness temperature fluctuations of the redshiftcd 
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21cm signal. The exact shape of the brightness tem- 
perature fluctuation power spectrum during reioniza- 
tion depends on the size and spatial distribution of the 
H II regions and should peak on the scale correspond- 
ing to the characteristic H II region size. On scales 
below this characteristic size the power in the fluctua- 
tions should decline. It is difficult to model this power 
spectrum without introducing considerable astrophysi- 
cal uncertainties, involving the formation and spectr a 
of the first stars and quasars l|Wvithe fc Loebl 12003) . 
feedback effect s that might impede th e formation of 
future sources l|Furlanetto fc Lo e^j2005 |). and shadow- 
ing by compact low mass halos (Barkan a fc Loeb1l2002l: 
Shap iro et alJl2004|) . Ultimately, hydrodynamical simu- 
lations with radiative transfer will be the most reliable 
method for following the evolution of the H II regions 
during the epoch of reionization. However, even these 
will inevitably be sensitive to uncertain details about the 
processes of star formation and quasar accretion which 
supply the ionizing photons. At present the large dy- 
namical range needed in order to follow reliably radia- 
tive transfer i n simulated boxes with periodic bound- 
ary conditions ( Barkana fc Loebll2004b() . is still orders of 
magnitude bey ond the values that state-of-the-art codes 
<IGnediril2004ISokasian et~all l2004) are able to achieve. 

In this paper we examine the polarization of the red- 
shifted 21cm radiation. The polarization could result 
from: (i) intrinsic properties of the sources; and (ii) pro- 
cesses that polarize the radiation as it travels towards 
the observer (secondary mechanisms). In the first cat- 
egory we will analyze processes which produce a non- 
isotropic population of the hyperfine triplet state, which 
upon spontaneous decay could p roduce 21cm polariza- 
tion. iCoorav fc Furlanettol l)2004) already considered the 
intrinsic polarization due to the Zeeman effect. They 
concluded that the effect was most likely too small to be 
observed for realistic values of the intergalactic magnetic 
fields at high redshifts. The second category of polariza- 
tion mechanisms is identical to the secondary processes 
which produce polarization in the Cosmic Microwave 
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Background (CMB) iZaldarriagal(l997l IHiJ 1200(1) . Here 
we will focus on the generation of polarization out of 
the free-streaming 21cm anisotropics through Thomson 
scattering by the reionized universe. 

The polarization anisotropics carry important new in- 
formation about reionization. For example, if the dom- 
inant source of polarization is Thomson scattering in 
the reionized universe, then the 21cm polarization would 
gauge the topology of ionized bubbles near the end of 
reionization as well as the optical depth to electron scat- 
tering after the bubbles overlaped. The latter mea- 
surement would help remove degeneracies among cos- 
mological parameters from the CMB data alone, and 
in particular substantiat e the inference abou t the opti- 
cal depth from WMAP ijSpergel et al.ll2003j) . The as- 
sociated constraints on the growth and morphology of 
H II regions during the epoc h of reionization ijZahn et al.l 
12003 IMcQuinn et aP fcOO.^ , can potentially remove the 
existing ten sion between the con clusions drawn from the 
CMB data iSpergel et al.ll2003D and the sizes of quasar 
H II r egions l|Wvithe fc Locb 2004a; Mcsinser & Haiman 
2004). 

An additional reason to study the polarization of 
the 21cm signal is that the expected low-frequency 
foreground is several orders o f magnitude br ighter 
than the redshifte d 21cm signal ijDi Matteo et alJl2002t 
IQh fc Mackl H003) . The foregrounds are expected to 
be smooth in frequency space, while the signal due to 
21cm emission should fluctuate between narrowly spaced 
bands, and so it was argued that the small signal can 
be extracted by taking the differences between maps at 
narrowly separated frequency bands ijZaldarriaea et all 
120041: IMorales fc Hewittl 12004}) . Since the foregrounds 
are expected to be polarized, a definite prediction for 
the polarization, both the signal and foreground, would 
further help us disentangle the signal from the fore- 
grounds. If the polarization signal is observable, it may 
be of a particular parity type (e.g. E-type or B-type, 
see iZaldarriaea fc Selialdll997l for a review of polariza- 
tion). This additional information could also be used to 
separate the signal from the foregrounds. Even if the 
polarization signal is not expected to be observable by 
a particular experiment, we can still use the prediction 
that any polarization detected by such an experiment 
should not correlate with the inferred signal in order to 
reduce the foreground contamination. 

The organization of the paper is as follows. In §2 we 
review the possible mechanisms which can produce a po- 
larized signal. In §3 we present numerical results for sec- 
ondary polarization anisotropies due to Thomson scat- 
tering in the reionized universe, which we expect to be 
the dominant source of polarization. Finally, §4 sum- 
marizes our main conclusions. We adopt the standard 
ACDM cosmo l ogica l model cosistent with the WMAP 
IjSpergel et all 12003]) data (0 6 = 0.044, fl m = 0.27, 
Ci v = 0.73, n = 1, erg = 0-9 and h = 0.72), and we 
show results for different values of the reionization opti- 
cal depth and the redshift of reionization. For simplicity, 
we ignore helium and assume a pure hydrogen plasma. 
We also focus on 21cm emission since the spin tempera- 
ture of hydrogen is likely to exceed the CMB temperature 
of Xcmb = 2.7 x (1 + z)K due to heating by X-rays; the 
heating is expected to occur soon after the first sources 
light up and long before the universe gets reionized, since 



it requires only < 10 2 eV (instead of > lOeV for ioniza- 
tion) per baryon. 

2. POLARIZATION MECHANISMS 

As already mentioned, the mechanisms that are ca- 
pable of producing polarization in the 21cm line divide 
into two broad classes: (i) intrinsic emission; and (ii) 
secondary mechanisms. In order to produce polarized 
intrinsic emission, the triplet state of the hyperfine split- 
ting needs to be excited in a non-isotropic way (e.g. 
through pumping by an anisotropic radiation field or 
splitting by a magnetic field which introduces a preferred 
direction) and then the spontaneous decay of this ex- 
cited state will produce polarized light. It is also nec- 
essary that collisional processes, which drive the system 
to thermal equilibrium, will not be more rapid than the 
spontaneous dec ay and anisotropic triplet pumping rates 
l|ZvgelmarJl2005|) . 

Even if the intrinsic emission is unpolarized, secondary 
mechanisms such as Thomson or resonant scattering 
along that linc-of-sight would produce polarization in an 
analogous manner to the secondary production of po- 
larization in the CMB. In this case, the polarization is 
sourced by an incident quadrupole brightness moment in 
the scatterer's frame of reference. The optical depth to 
electron scattering i nferred by WMAP r ~ 0.17±0.04 
ijSpergel et al.ll2003|) . implies that secondary scattering 
could potentially produce a non-negligible effect. The 
fundamental difference between the scattered 21cm po- 
larization and reionization bump for the CMB polariza- 
tion is the source of the incident quadrupole. The in- 
cident quadrupole in the CMB case originates at the 
surface of last scattering at the redshift of recombina- 
tion (z — 10 3 ) and is due to the projection and free- 
streaming of the correlated CMB temperature fluctua- 
tions. The free-streaming distance is always the horizon 
size as seen by the scatterer. While the basic physical 
process is identical, these details are changed in the 21 cm 
case. First, the free-streaming distance from the emitter 
to the scatterer is not always the horizon size, but set by 
the frequency of the observed line radiation. The dom- 
inant sources of the incident quadrupole during reion- 
ization are likely to be baryonic density fluctuations in 
the H I gas and ionization fraction fluctuations based on 
the topology of H II regions. The optical depth for reso- 
nant 21cm absorption is typically smaller th an a percent 
(|Carilli et alJl200l IPlIrlanetto fc Loebll2002|) and will be 
ignored in our discussion. 

2.1. Intrinsic Emission 

Isotropic symmetry needs to be broken in order pro- 
duce net polarization (i.e. a preferred direction) in 
21cm line emission from an ensemble of hydrogen atoms. 
There are two simple ways to break isotropy: (i) Zee- 
man splitting of the triplet state due to a directional 
magnetic field; and (ii) Wouthuysen-Field effect due 
to an anisotropic Lya radiation field. The possibility 
of the Zeeman effect has already been considered by 
ICoorav fc Furlanettol ((2004) . In order for the signal to be 
measurable by SKA, a coherent field strength of > lOO^G 
on Mpc scales is required at z ~ 6-10 , with an energy 
density comparable to the CMB at that time. Such a 
fiel d strength is well above current Faraday-rotation lim- 
its ifVallcc 2004) on the intergalactic field outside galaxy 
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clusters (keeping in mind the adiabatic enhancement fac- 
tor of (1 + z) 2 to high redshifts), and is orders of mag- 
nitude a bove theoretical expectat i ons from astrophysic al 
sources l|Furlanetto fc Loefjll200ll iKulsrud et alJ l!997). 

The Wouthuysen-Field effect is the mechanism which 
couples the gas temperature to the hyperfine spin tem- 
perature through the absorption and e mission of Lya 
photons (| Wouthuvsenl Il95l iFieldl Il959|) . There is a fi- 
nite probability for a 2p state, which was reached by Lya 
absorption by the singlet Is state, to spontaneously de- 
cay into the triplet Is state. This process will convert a 
fraction of the singlet states to triplet states. When the 
Lya scattering rate is high, the spin temperature which 
determines the relative populations of the singlet and 
triplet states, will equal the kinetic temperature of the 
gas (which sets the "color temperature" of the Lya pho- 
tons) . During reionization the kinetic temperature of the 
gas is heated (e.g. by a small level of X-rays produced by 
the first compact objects) to a value much greater than 
the CMB temperature. The corresponding increase in 
the spin temperature allows the 21cm line to be seen in 
emission during the late stage of reionization. 

If the Lya radiation as seen by a hydrogen atom has a 
quadrupole moment, then the 21cm line emitted by this 
atom will be pol arized (as long as atomic collisions can 
be neglected; see IZvgelmanll2005j) . The basic idea is that 
continuum radiation from the first high-redshift sources 
will be rcdshifted into the blue wing of the Lya line pro- 
file. This radiation field could potentially be anisotropic 
because the ionizing sources ar e highly clustered at high- 
redshift ijWvithe fc Loebl2005j) . Unfortunately, the large 
optical depth near the Lya resonance (~ 10 5 for a neu- 
tral IGM) quickly isotropizes the incident radiation as it 
is redshifted through the line profile by the cosmological 
expansion. Therefore only the photons in the far blue 
wing of the Lya line profile are expected to produce a 
non-isotropic pumping of the triplet state distribution. 
However, these photons are much less likely to scatter 
than the isotropized photons near the center of the line. 
Even in the vicinity of a bright quasar, the anisotropy of 
the incident radiation would induce polarization only for 
the neutral IGM within a thin skin outside the quasar's 
H II region whose spatial width is of order the mean- 
free-path at t he center of the thermal ly-broadenned Lya 
line, < 0.1 pc l)Wvithe fc Loebl2004c]) . We conclude that 
pumping by an anisotropic Lya radiation field is unlikely 
to produce observable intrinsically polarized 21cm line 
radiation. 

2.2. Secondary Mechanisms 

Thomson scattering by free electrons along the line- 
of-sight can polarize intrinsically-unpolarized emission if 
the radiation field as seen by the scattering electrons 
possesses a quadrupole moment. The inhomogeneous 
distribution of the emitting H I gas naturally gener- 
ates a 21cm quadrupole as the radiation free-streams to- 
wards the free electrons in the reionized universe. The 
21cm emission originates from the neutral fraction of 
the IGM, and so fluctuations in the ionization fraction 
will also produce fluctuations in the 21cm brightness. 
There should be fluctuations in the redshift of reion- 
ization along d ifferent directions due to large scale in- 
homog eneities llBarkana fe Loebl l2004at |Wvithe fc Loebl 
l2004cD . Unfortunately, the detailed correlation function 



and size distribution of H II regions depends strongly on 
uncertain details of the astrophysics of star and black 
hole formation. This dependence introduces inherent 
uncertainties into the theoretica l modelling of the ion- 
ization fraction powe r spectrum ijFurlanetto et alJl2004t 
McQuinn et al.ll2005j) . 

The end of reionization, however, offers a special 
regime where robust model-independent predictions can 
be made. This was demonstrated by Wyithe & Loeb 
(2004b) who derived the characteristic size of individ- 
ual H II regions just before their final overlap based on 
the extended Press-Schetcher theory. This derivation, 
which was based on the general physical considerations of 
cosmic variance and light propagation delay, eliminates 
much of the dependence on the astrophysical details. We 
will therefore incorporate this more robust approach and 
derive the 21cm polarization power spectrum from it. For 
simplicity, we will assume that Thomson scattering starts 
after the final overlap of the H II regions, since most of 
the Thomson scattering optical depth is naturally in the 
regime where the filling fraction of ionized regions is sub- 
stantial. In our model, the 21cm quadrupole is sourced 
by the Poisson fluctuations of the H II regions at the 
end of reionization. This effect will produce the largest 
signal, however for completeness in §3 we will also dis- 
play results from the baryon density fluctuation power 
spectrum. 

The two Stokes' parameters Q and U of polarized radi- 
ation can be decomposed i nto the rotationally invarian t 
E-type and B-type fields (jSeliak fe Zaldarriagai ri996). 
Due to parity invariance, perturbations which rotation- 
ally transform as scalars and therefore lack any handed- 
ness, will only produce E-type polarization. The 21cm 
quadrupole is sourced by baryon density and ionization 
fraction fluctuations and should only produce E-type po- 
larization since the induced polarization is azimuthally 
symmetric about the initial perturbation wavevector. 
This is an important conclusion because it can be used to 
eliminate polarized foreground emission and instrumen- 
tal systematics. 

Pr opagation effects, such as gravitational lens- 
ing IZ aldarriaga fc Seliakl 119981) and Faraday ro tation 
l|Kosowskv fc Loebl 1199a iKosowskv et alJ I2005j) . will 
convert the E-type polarization produced by Thomson 
scattering into B-type polarization. The amplitude of 
B-modes prod uced by gravitational lensi ng should not 
be significant ijZaldarriaga fc Seliakl 1199^1 . however the 
effect of Faraday rotation may be quite substantial at 
the relevant low frequencies. It is difficult to estimate 
the properties of the appropriate rotation measures since 
they depend on the magnetic field amplitude and co- 
herence structure, not just in the Galactic interstellar 
medium, but also in the high redshift IGM. Nevertheless 
the effect of Faraday rotation has a known frequency de- 
pendence and it is possible that the developed techniques 
for foreground removal may be adapted to address this 
challenge. 

While calculating the effect of Thomson scattering on 
the incident 21cm quadrupole we will ignore the inho- 
mogeneites present in the distribution of scattering elec- 
trons. For simplicity, we will calculate the expected po- 
larization signal assuming a sudden uniform reionization. 
In reality, prior to the end of reionization the ionization 
fraction is less than unity, with full ionization only in lo- 



4 



calized H II regions around biased clusters of sources. 
Since the ionized bubbles contain free electrons, they 
would also scatter the incident quadrupole and produce 
polarization. The resulting visibility function will be in- 
homogeneous and scattering will be able to produce B- 
mode polarization. In all calculations of the effect for 
the CMB, the produced B-mode power spectrum was 
confined to small scales and was several orders of mag- 
nitude small er in amp litude than the large scale E-mode 
polarization (Hu 2000); we do not expect this conclusion 
to differ for the 21cm line radiation, and so our working 
assumption is that the large scale polarization only de- 
pends on the optical depth of the homogeneous universe. 

In §3 we will demonstrate that the effect of assuming an 
instantaneous, instead of gradual, homogeneous reioniza- 
tion history is modest by showing that the polarization 
power spectrum can be well approximated by a polariza- 
tion power spectrum calculated from an instantaneous 
reionization history. Of course, the two power spectra 
are not identical; however, the appropriate power spec- 
trum error bars, both cosmic variance and instrumental 
noise, must be included in making the comparison. Si- 
miliarly it has been noted that the reionization history 
can be approximated as instantaneous for low sensitivity 
CMB polarization measurements (such as WMAP), but 
for upcoming nearly cosmic variance limited measure- 
ments (suc h as Planck) a better model of reionization i s 
necessary ijHaiman fc Holderll2003t IHolder et alJ 12003). 
Since the introduced theoretical uncertainties should be 
modest compared to the large instrumental noise and 
foreground contamination, we will use the assumption of 
instantaneous reionization in this paper. 

As customary, we write down the polarization 

transfer function i n th e line-of-sight fo rmalism 

ijSeliak fc Zaldarriagalll99(l IHu fc White! Il997l) . which 
is the Legendre expansion of the E-type polarization 
induced by a single Fourier mode perturbation, 



where k is the primordial wavevector, / is the observed 
Legendre mode, and rj is conformal lookback time, de- 
fined as 

Here we also define the visibility function, which is the 
probability that a photon last scattered at "q, as 

9(V) = -r^ r(v \ 0) 



drj 



where 



dr 
drj 



a T x e (z)n b (l + z) 2 , 



(4) 



<jt is the cross-section for Thomson scattering, x e (z) is 
the ionization fraction and is the average baryon num- 
ber density today. The incident quadrupole on the scat- 
terer is n(fc, rj) and v is the observed frequency of the 
21cm line. Since Thomson scattering is achromatic, the 
change in the observed frequency will simply be due to 
the cosmological redshift, thus we can relate the observed 
frequency to the rest frame frequency, vq — 1.4 GHz, as 
v = i^o/(l + z e), where ze is the emission redshift. 



The smooth intergalactic H I gas has a 21cm 
optical depth fFleldl IT959t iBharadwai fc Ali I l200l 



iBarkana fc Loebll2005a|) 

_ 3c 3 hA 10 x H n b (l + S)(l + zf 
T ~ Wk B vlT s H{z) 



(1 + z) dv r 
H(z) ~dr 



(5) 

where A w is the spontaneous emission coefficient of the 
hyperfine transition, Tg is the spin temperature that de- 
termines the relative populations of the triplet and sin- 
glet states, xh is the neutral fraction of the hydrogen 
gas, and S is the baryonic overdensity. Here dv r /dr is 
the radial gradient of the linc-of-sight peculiar velocity 
of the H I gas. The optical depth is proportional to the 
path-length over which the a 21cm photon resonates with 
hydrogen atoms, which is in turn inversely proprtional 
to the rate by which the Doppler shift of the medium 
changes along its path. Equation contains contribu- 
tions to the optical depth from both the cosmological 
expansion and the peculiar velocity of the gas. The fluc- 
tuation in the brightness temperature, calculated from 
the equation of radiative transfer, is 



™ , , Ts - Tomb 

blb{z) = — ; - — T, 



(1 + *) 

23[l + (l + ti 2 )5}x H 



10 



1/2 



x {Ts -J CMB) ( 1±£ 1 mK, 



(0 



(7) 



where fi is defined as /x = k ■ n, the cosine of the angle 
between the wavevector of the primordial perturbation 
and the line of sight. Since the secondary polarization 
mechanism is achromatic, we will simply consider the di- 
mensionless brightness temperature fluctuation, ip, which 
is defined as 



ip = [1 + (1 + ^)5]xh 



(Ts — Tomb) 
Ts 



(8) 



In the limit Ts 3> Tomb, which we expect to hold late in 
reionization, this can simply be expressed as ip = xh [1 + 
(1 + ^)5}. 

These correlated brightness temperature fluctuations 
will produce a quadrupole incident upon the scatterers 
because free-streaming of fluctuations tr ansfers power 
from low Legendre modes to higher ones l)Hu fc White 
1997). Performing this projection and free-streaming, 
the incident dimensionless qua drupole can be expressed 
as ()ZaldarriagallT997l 151120001 



Tl(k,r),is) 



+ i 



(9) 



= .h[k(VE - T])]x H (k,i] E )(l + S(k,j] E )) 
- sSWve - r))]x H (r) E )5(k,r) E ), (10) 

where r] E is simply the conformal distance to redshift 
of emission, z E . Each prime denotes a derivative of the 
spherical Besscl function; the derivative term in equation 
(10) is produced by the peculiar velocity term in equation 
©■ 

In difference from the calculation of CMB anisotropies, 
we need to define two separate transfer functions for the 
polarization produced by the baryonic density fluctua- 
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tions, 

AE „ . 3 1(1 + 2)1 f m , g{n) . s 

' A f(*.-)=4vMry 

xljat^-TjJl-i^B-tj)]], (11) 
and the ionization fraction fluctuations, 



We must make this distinction since we cannot simply 
model the ionization fraction fluctuations in terms of the 
baryonic density fluctuation or primordial curvature per- 
turbation. 

The observed E-type polarization power spectrum can 
be expressed as 

c ?(y) = ~ / k " dk tfkPsW [ s Af(k,u)} 2 

+P x (k)[ x A?(k,v)}% (13) 

where Pg(k) is the baryonic power spectrum and P x (k) 
is the ionization fraction power spectrum. Here we have 
defined xr to be the average ionization fraction. Wc 
will ignore the cross-correlation between 5 and xh ■ The 
baryon density fluctuation can be expressed as 

5(k, VE ) = T(k)D( VE )((k), (14) 

where T(k) is the standard transfer function calculated 
by CMBFAST 8 , D(rf) is the linear theory growth func- 
tion and C(fc) is the primordial curvature fluctuation pro- 
duced by inflation. Thus, the baryonic power spectrum 
can be expressed as 

P 5 (k,r) E ) = T 2 (k)D 2 (r] E )P c (k), (15) 

where P$(k) is the scale-invariant primordial curvature 
power spectrum. 

As previously mentioned, there will also be correla- 
tions in the ionization fraction fluctuations. The rela- 
tive importance of the contributions from the baryonic 
power spectrum and the ionization fraction power spec- 
trum depends on the mean ionization fraction as well 
as the topology of H II regions. The calculation of 
P x {k) is complicated by uncertainties about the ionizing 
source population and the radiative transfer of the ioniz- 
ing radiation. There have been recent attempts to semi- 
analyti cally model this power spectrum based on the halo 
model iJMcQuinn et al.ll2005j) . which divides the power 
spectrum into contributions from the internal distribu- 
tion of individual halos and the large scale clustering of 
halos that are biased tracers of the linear ba ryonic mat- 
ter power spectrum l|Coorav fe Shethl |2002fl . The con- 
tribution from individual halos introduces Poisson fluc- 
tuations in the H II regions and dominates the power 
spectrum on small scales. On large scales, the halos are 
a biased tracer of the linear baryonic power spectrum. 
The transition between these two regimes occurs near 
the characteristic bubble size. 

In this paper we focus on the contribution to the inci- 
dent quadrupole from Poisson fluctuations in the number 

8 http://www.cmbfast.org/ 
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Fig. 1. — A schematic illustration of the scattering geometry. 
Each electron along the dotted (rightmost) line-of-sight scatters 
21cm radiation that was emitted just before the end of reionization 
at a redshift zr from a spherical shell around it. All emitting 
shells are tangent to the bounding shell around the observer. 21cm 
radiation emitted by the patchy H I distribution in these shells and 
scattered after reionization by the electrons at the shell centers will 
show polarization at a wavelength of 21cm X (1 + Zr). The two 
solid lincs-of-sight on the left which are separated by an angle 9 will 
show correlations in their polarization amplitudes due to scattering 
electrons with intersecting emission shells. 



of H II regions. For simplicity and definitiveness we will 
only include this effect at the end of reionization, namely 
the surface of bubble overlap, where we can robustly cal- 
culate the size and number density of H II regions based 
on the general cosid erations of cosmic varia nce and light 
propagation delay 9 ijWvithe fe Loebll2004b|) . As long as 
the volume filling fraction of ionized bubbles is small, the 
21cm fluctuations are small, and so most of the signal 
naturally originates around the time of bubble overlap 
anyway. We will ignore the correlations of the ionization 
fraction due to large scale biasing of the H II regions, 
as well as, cross correlations between the baryon den- 
sity and the ionization fraction. Biasing effects can only 
be important on scales above the characteristic size of 
the H II regions. However at the end of reionization the 
characteristic size of t he H II regions is extre mely large 
(~ 70 comoving Mpc) ijWvithe fc Loebl2004b]) . On these 
scales the amplitude of the baryonic power spectrum is 
quite small so the biased H II region power spectrum 
will be modest at the end of reionization. Besides these 

9 Note that inclusion of the light propagation delay is essential 
for the calculation of the size of an H II bubble at the end of 
reionization, as viewed from the location of a scattering electron. 
Any calculation of the typical bubble size in a spacelike snapshot of 
the universe (as commonly done in the literature) is not adequate 
for this purpose. 
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effects are primarily important on large angular scales, 
where the conversion of the temperature brightness fluc- 
tuations to the necessary quadrupole anisotropies is sup- 
pressed; therefore, we will find them to be below realistic 
instrument detection thresholds. 

The observed fluctuations in the ionization fraction are 
caused by fluctuations in the local number of H II regions. 
Poisson statistics of point-like objects is described by the 
spatial two-point correlation function, 



(5x(ri)5x(r 2 )) 



1 



n(ri) 



£ (3) (n-r 2 ), 



(16) 



where n is the average number density of H II regions 
and Sx is the fluctuation in the ionization fraction. Since 
the H II regions have a finite size, we must convolve their 
density field with a window function 



Sx(ri 



= J d^r 2 W R {\rx -r 2 \)6x(r 2 ), 



(17) 



where we assume a Gaussian window function for com- 
putational simplicity 



W R (r) 



The two-point correlation function is given by 
(Sx(ki)5x*(k 2 )) = 



(18) 



3... W 3 ^. a iki-ri _—• Sfei-n / j3„, j3 



<fr 1 d i r 2 e 



drvuryi 



\r 1 -y 1 \ 2 /2R 2 - \r 2 -y 2 \ 2 /2R 2 



Simplifying we find 



-{5x{y x )8x{y 2 )tlQ) 



(6x(k 1 )Sx*{k 2 )) = 

4 / d 3 ri d 3 r 2 e lki ri < 
n J 



-ikx -T"i 



-D 2 /B 2 



(20) 



where we have defined D = (ri — r 2 )/2. As R — > 
0, the Gaussian window function approaches (5^ 3 '(ri — 
r 2 ) and then we find P x (k) — 1/n, where we define 
the ionization fraction power spectrum, P x (k), through 
{5x(kt)Sx*(k 2 )) = (27r) 3 5( 3 )(fci - k 2 )P x (k 1 ). We can 
now analytically calculate the two-point function for fi- 
nite bubble size 



-k A R 2 



(21) 



The translational invariance of (8x(ri)Sx(r 2 )) guar- 
antees that (Sx(k%)Sx* (k 2 )) will be proportional to 
(5' 3 '(fci — k 2 ). The cutoff in the temperature bright- 
ness fluctuation power spectrum due to the finite bubble 
size will cause a corresponding decay in the polarization 
power spectrum for multipole index I > r/ji/R. Below 
this scale, the afforementioned signal due to baryon fluc- 
tuations will come to dominate. However, as we will see 
in the next section, this contribution to the signal is far 
below the instrumental noise detection threshold. There- 
fore we are justified in ignoring it. 

We assumed n(r) to be constant because of the large- 
scale homogeneity of the universe; the coordinate r is the 
comoving distance on the hypersurface of simultaneity 



corresponding to the time of reionization, not a coordi- 
nate labeling our past lightcone. This Poisson fluctuation 
contribution to the polarization power spectrum corre- 
sponds to the intersection (within the region defined by 
the bubble window function) of the past light cones of the 
two electron scatterers, which both lie on our past light 
cone, on the hypersurface of simultaneity corresponding 
to the end of reionization. Figure [I] schematically shows 
this effect. On the righthand side of the figure, scatter- 
ers along a given line-of-sight see 21cm radiation from 
progressively larger shells because the elapsed time be- 
tween emission and scattering for an electron increases 
for electrons closer to the observer. On the lefthand side 
of the figure, the shells of 21 cm radiation emisson are 
displayed for two electrons along two lines-of-sight sepa- 
rated by an angle 9. These shells are distinct from the ob- 
server's shell because they correspond to past light cones 
of the electrons which intersect the hypersurface of si- 
multaneity when reionization ended at different physical 
positions. In the figure, the shells for the two electrons 
(red and blue) intersect and therefore the polarization 
produced along the two lines-of-sight will be correlated 
due to Poisson fluctuations. 

At the end of reionizatio n, which we define to b e the 
surface of bubble overlap ijWvithe fc Loebll2004bj) . the 
mean distance between H II regions is equivalent to the 
characteristic size of the the bubbles. Thus n, the av- 
erage bubble density, equals 1/V, where V = 47ri? 3 /3 
is the characteristic bubble volume. This is only true 
at the end of reionization; earlier on, during the process 
of reionization, the characteristic bubble size is smaller 
than the mean bubble separation. For our model, we find 
that the polarization power spectrum is given by 



o n3 

Cf(y) = — I k' z dk 



-R 2 k 2 



UAf(fc^) 



(22) 



Since the brightness temperature fluctuations are the 
source of the quadrupole which produces polarization 
through Thomson scattering, the polarization and tem- 
perature signals should be correlated. Actually only the 
E-type polarization is correlated with the temperature 
fluctuations, as parity prevents correlati ons with B-type 
polarization (jZaldarriaga fc Selialdll997j) . Depending on 
the characteristics of reionization, this cross-correlation 
signal may be much easier to detect. For this reason the 
CMB temperature-polarization cros s-correlation spec- 
trum was firs t published by DA SI ijKovac et alJ 12002^) 
and WMAP ijKoeut et aQ Eoog). It is straightforward 
to adapt the above formalism. The temperature trans- 
fer function is again divided into contributions from the 
baryon density fluctuations 



(23) 



<yAf(fc, v) = ji(kr} E ) - j"(kr) E ), 

and the ionization fraction fluctuations 

x Aj(k,v)=j l (kr lE ). (24) 

Then the observed temperature polarization cross- 
correlation spectrum can be written as 

C ^ v ) = \\ k2dk lx 2 H Ps(k)sAF(k,is) s Aj(k,v) 

+P x {k) x Af(k,p) x Aj(k,u)}. (25) 

We have ignored effects due to the spectral response 
functions of the observing instruments. These are 
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Fig. 2. — The brightness temperature (red, dashed) and polariza- 
tion (black, solid) due to baryon density fluctuations for emission 
redshift of 2£ = 20 and reionization redshift zr = 17. 
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Fig. 3. — The polarization power spectrum sourced by baryon 
density fluctuations for an emission redshift ze = 20 and reioniza- 
tion redshift zr = 17 (black, solid), zr = 15 (red, dotted), zr = 10 
(blue, dashed) and zr = 6 (green, long dashed). 



straightforward to incorporate into our formalism. Since 
we will find that the dominant signal, which is produced 
by H II region Poisson fluctuations, is rather feature- 
less and exponentially damped below scales of tens of 
arcminutes (/ > 10 3 ) reasonable instrument bandwidths 
(Af ~ 0.4 MHz) will not substantially affect our results. 

3. RESULTS 

Based on the formalism developed in §2, we present 
numerical results for a variety of emission and reioniza- 
tion redshifts in order to understand how the signal de- 
pends on these parameters. We consider the simplest 
model in which reionization is assumed to be sudden, 
so the ionization fraction is a step function, and homo- 
geneous, except for the effect of Poisson fluctuations of 
H II regions on the incident 21cm quadrupole. This im- 
plies that the 21cm power spectrum will simply be due 
to correlations in the baryonic density and that the gas is 
always completely neutral, xh = 1, when the redshift of 
the emission is greater than the redshift of reionization. 
When we observe emitted radiation from the redshift of 
reionization, we include the contribution from the Pois- 
son fluctuations in the H II region distribution. Once the 
universe becomes reionized, it does so uniformly. Our ap- 
proximate treatment ignores scattering due to the patch- 
iness of the universe while it is only partially ionized. 
We adopt for our first treatment this simplified model 
because it has the smallest number of free parameters; 
more complicated models with uncertain astrophysical 
parameters can be considered in the future. 

Figure [5] compares the temperature brightness power 
spectrum (Cf) and the polarization power spectrum 
(Cf) due to baryon density fluctuations for ze = 20, 
Zr = 17. The features in Cj and Cf are quite different. 
In addition to being at a lower amplitude, Cf peaks and 
begins to oscillate and decay at high I. The polarization 
transfer function, defined in equation JQ, is the line-of- 
slight integral of two spherical Bessel functions, weighted 
by the visibility function and additional factors due to 



the polarization spin lowering operators. The properties 

of the spherical Bessel function are well-known: when 

x <C 1 ji (x) ~ x , it has a peak at x ~ I and then at large 

x it oscillates and decays as ji(x) ~ sin(a; — irl/2)/x. We 

will utilize these features to understand the behaivor of 
r<E 

Figure y] shows the polarization power spectrum for 
various reionization redshifts, Zr = 17,15,10,6 (top-to- 
bottom), and a single emission redshift, ze = 20. The 
basic physical processes become clear by comparing these 
various power spectra. We will first consider the angular 
position of the peak in the polarization power spectra. 
The spherical Bessel function ji(krj) peaks at I ~ kr\ 
and likewise j2 [k(f]E ~ v)] a ^ 2 ~ k(rjE — rj). These two 
constraints imply that the polarization power spectrum 
peaks at 

Ipeak ^ ) (26) 

VE ~ VR 

where r\E and t]r are the comoving distances to the red- 
shifts of emission and reionization, respectively. Here 
we have effectively assumed an infinitely sharp visibil- 
ity function. The inclusion of the altered free-streaming 
term, due to peculiar velocities, only has a minor effect 
on this argument. Thus changing zr, for a given ze, 
will shift the power spectrum features to higher I. This 
is clearly observed in Fig. and numerically the peak 
locations agree with equation (|2l)|) . 

The oscillatory features are due to the free-streaming of 
the monopole brightness temperature fluctuations. The 
incident quadrupole on a scatterer, equation JjJJ, con- 
tains the spherical Bessel function, j2[k(r/E — T))], which 
describes how a monopole fluctuation of wavevector k at 
t\e becomes a quadrupole fluctuation at r\. The pro- 
jection of this oscillatory function, the free-streaming 
quadrupole, on the sky causes the features in Fig. 0. 
These are not the analogs of the acoustic peaks observed 
in the primary tempera ture anistotropies of the CMB 
ijBarkana fc Loeb1 l2005c) . Both are caused by projecting 



an oscillatory function on the sky; however the acoustic 
oscillations originate from pressure waves in the baryon- 
photon fluid prior to recombination. 

The decay of Cf at high I is explained by oscilla- 
tory cancellation of the line-of-sight integration. The 
two spherical Bessel functions would have had different 
phases and therefore integrated to zero if the range of in- 
tegration was infinite. When the wavelength of the per- 
turbation becomes comparable in size to the width of the 
visibility function, this effect begins to become impor- 
tant. Since the visibility function becomes progressively 
more peaked at higher reionization redshift, progressively 
smaller intervals along the line-of-sight contribute to the 
transfer function and therefore there is less damping. We 
observe this phenomenon in Fig. |3 as the high I tail de- 
cays less and retains oscillatory features for the curves 
corresponding to the higher reionization redshifts. If the 
change in the ionization fraction is not sudden as we have 
assumed but rather gradual, then these oscillatory fea- 
tures would be smoothed out. 

In these examples we have assumed that reionization is 
inst antaneous, while in reality reionization must be grad- 
ual l|Furlanetto fc Loebl \mm . As mentioned in §2.2, 
the adopted reionization history will imprint a signature 
on the polarization power spectrum. Next we demon- 
strate that the polarization power spectrum sourced by 
baryon density fluctuations and produced for a given 
gradual reionization history can be well approximated 
(within the appropriate cosmic variance and instrumen- 
tal noise errorbars) by an instantaneous reionization his- 
tory (with the total optical depth not being necessarily 
the same in the two reionization histories). In the bottom 
panel of Fig. 0] the polarization power spectra sourced 
by baryon density fluctuations are displayed for the in- 
stantaneous and gradual reionization histories plotted in 
the top panel. While there are differences between the 
two reionization histories, they are fairly modest com- 
pared to the much greater instrumental noise and model 
uncertainties. The use of the simple model of instanta- 
neous reionization is therefore well justified for the scope 
of this paper. 

Figure [S] displays the polarization power spectra 
sourced by Poisson fluctuations of H II region, for differ- 
ent redshifts of emission and reionization: ze — zr = 30 
(black, solid), ze = zr = 25 (red, dotted), ze = zr = 17 
(blue, dashed), ze = Zr = 10 (black, long dashed) and 
ze = zr = 6 (red, dot dashed). The theoretical detection 
threshold of SKA is also shown for a 1 month integration 
time and a 1 year integration time. We can understand 
the basic features of the polarization power spectra by 
inspecting equation (|22|) . The amplitude of the polariza- 
tion power spectra is proportional to the characteristic 
volume of an H II bubble at the time of bubble overlap. 
If reionization completes at a higher redshift, then the 
H II bubbles are smaller and therefore more numerous 
at the surface o n bubble overlap which defines the end 
of reionization l)Wvithe fc Loebl I2004H) . Subsequently, 
this reduces the level of the Poisson fluctuations, which 
is inversely proportional to the H II region number den- 
sity. The window function related to the finite size of 
a typical H II region, produces an exponential decay in 
the power spectrum on smaller spatial scales. Thus the 
power spectra corresponding to higher reionization red- 
shifts, and therefore smaller characteristic sizes, begin to 
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Fig. 4. — Top Panel: The gradual and instantaneous reionization 
histories appropriate for the polarization power spectra displayed 
in the botton panel. Bottom Panel: The polarization power spectra 
are displayed for the instantaneous and gradual reionization histo- 
ries assuming that the 21cm fluctuations are sourced by density 
inhomogeneities. The color coding is the same in both panels. 
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Fig. 5. — The polarization power spectrum sourced by Poisson 
fluctuations of H II regions for an emission and reionization red- 
shift = zr = 30 (black, solid), ze = 2 R = 25 (red, dotted), 
Z E = Z R = 17 (blue, dashed), = zr : = 10 (black, long dashed) 
and ze = zr = 6 (red, dot dashed). The theoretical detection 
threshold of SKA is also shown for a 1 month integration time and 
a 1 year integration time with a bandwidth of 0.4 MHz. 



exponentially decay at higher values of I. 

In Fig. [SJthe amplitude of the polarization power spec- 
tra increases with increasing reionization redshift. This 
conflicts with the naive expectation that the polariza- 
tion power spectra will be smaller at high redshift since 
the Poisson fluctuations are smaller. The reason is that 
the decrease in the characteristic H II region size is rather 
gradual and the visibility function is significantly increas- 
ing at high redshifts. This is the main reason for the fact 
that the polarization power spectra with low reionization 
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Fig. 6. — The temperature-polarization cross-correlation spectra 
sourced by Poisson fluctuations of H II regions for an emission and 
reionization redshift = zr = 30 (black, solid), zg = zr = 25 
(blue, long dashed), zg = zr = 17 (red, dashed), z^ = zr = 
10 (black, dotted) and ze = zr = 6 (red, dot dashed). The 
theoretical detection threshold of various instruments : SKA - 1 
month, LOFAR - 1 year and MWA - 1 year of integration time all 
with a bandwidth of 0.4 MHz. 



rcdshifts have lower amplitudes. At yet higher redshifts 
a turnover does occur and the decrease in characteristic 
H II region size overcomes the increase in optical depth. 

The polarization power spectrum detection thresh- 
old due to instrument noise is included on Fig. [3] 
for the Square Kilometer Array (SKA), which is the 
next generation low frequency radio interferometer. We 
adopt the power spectru m detection threshold model of 
iZaldarriaea et al.1 l)2004|) with updated parameters. We 
assume an instrument bandwidth of Av ~ 0.4 MHz 
throughout the paper; increasing the bandwidth will de- 
crease the instrument detection threshold as At^ 1 / 2 . In 
reality, the relevant power spectrum error bars for a ra- 
dio interferometer are a complcated function of both the 
exact instrumental design and analysis method, so the 
detection thresholds in this work should be viewed as a 
crude estimate (assuming that the foregrounds can be 
properly removed). Measurements of the proximity ef- 
fect around high redshift quasars implies that the end of 
reioniz ation, i.e. the surface of bubble overlap, is near 
z = 6 l|Wvithe & Loebll2004a|) . We find that the polar- 
ization signal from even a redshift as low as ze — Zr = 6 
is detectable with a 1-year integration time. A higher 
redshift of reionization will only increase the amplitude 
of the polarization power spectrum. 

In §2 we predicted that the temperature polarization 
cross-correlation spectrum should be significantly larger 
than the polarization power spectrum. Since the polar- 
ization power spectrum may only be detected by a year 
long integration with a next generation instrument, we 
will compute the cross-correlation spectrum and compare 
its amplitude to the noise detection thresholds of up- 
coming instruments. Figure El displays the temperature- 
polarization cross-correlation spectra, sourced by H II re- 
gion Poisson fluctuations, for different redshifts of emis- 



sion and reionization: ze — Zr — 30 (black, solid), 
ze = zr — 25 (blue, long dashed), ze = zr — 17 (red, 
dashed), ze — zr — 10 (black, dotted) and ze — Zr — 6 
(red, dot dashed). Also shown are the instrument de- 
tection thresholds for SKA - 1 month integration time, 
LOFAR - 1 year and MWA - 1 year. We can understand 
how the features of these spectra differ from the polariza- 
tion power spectra of Fig. |5]by comparing the formula for 
the temperature polarization cross-correlation spectrum 
(Eq. I25|) with the formula for the polarization power 
spectrum (Eq. l2*2*j) . 

There are two obvious differences between the spec- 
tra. The temperature polarization cross-correlation spec- 
tra have larger amplitudes and shallower slopes on large 
scales. Clearly the polarization transfer function (Eq. 
^) has a smaller amplitude than the temperature trans- 
fer function (Eq. |23J because of the low optical depth, 
as well as, the effects of the quadrupole free-streaming 
and conversion to E-type polarization. The steep large 
scale slope observed in the polarization power spectra 
(Figs. EEt results from the quadrupole nature of the 
production of polar ization through Thomson scattering 
ijHu &: White! 11 997). Each polarization transfer func- 
tion roughly scales as I 2 on large scales, so the polar- 
ization power spectrum should naturally be steeper than 
the temperature polarization cross-correlation on these 
scales. Once again, on small scales we observe the same 
exponential decay caused by the H II region window func- 
tion. We find that even if reionization ended at Zr = 6, 
this signal should be easily detectable by SKA and very 
close to the detection threshold for LOFAR and MWA. 
If new observations of high redshift quasars imply that 
the surface of bubble overlap occured at a higher red- 
shift, then the detectability of the temperature polar- 
ization cross correlation signal increases for LOFAR and 
MWA. The result is very encouraging because it is based 
on simple analytical considerations which are robust to 
variations in astrophysical parameters. These uncertain 
parameters are all lumped into the value of zr, in our 
formulation. 

4. DISCUSSION 

Figure [S] implies that the future Square Kilometer 
Array (SKA) will have sufficient instrument sensitivity 
to detect the temperature polarization cross-correlation 
power-spectrum of redshifted 21cm fluctuations. The 
forthcoming Mileura Widefield Array (MWA) and the 
Low Frequency Array (LOFAR) might also have sensitiv- 
ity to detect the predicted signal. If new observations of 
high redshift quasars push this redshift of bubble over- 
lap to higher values then the possibility of a detection 
increases. The practical feasibility of such a detection 
will be better known within a few years, as soon as the 
first polarized foreground maps will be produced by LO- 
FAR or MWA. The instrumental detection thresholds 
used here are based on simplifications and should only 
be interpreted as rough estimates. However, the results 
shown in Figure portray the optimistic forecast that 
the polarization signal is sufficiently large for it to be 
detectable. 

We identified the dominant producer of polarization 
to be Thomson scattering of an incident quadrupole 
moment of the 21cm radiation. The 21cm quadrupole 
is sourced by the free-streaming of brightness fluctu- 
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ations from either correlated baryonic density fluctua- 
tions or Poisson fluctuations of H II regions. The ampli- 
tude of these Poisson fluctuations can be robustly pre- 
dicted based on the general considerations of cosmic vari- 
ance and light propagation delay at the end of reion- 
ization. We have found that for a reionization redshift 
of zr — 6, which is suggested by the observations of 
the Gunn-Peterson effect, the temperature polarization 
cross-correlation spectrum amplitude due to Poisson fluc- 
tuations should be of order ~ 0.1 mK and the polariza- 
tion power spectrum amplitude ~ 3 fj,K. Strong cross- 
correlation between the brightness temperature and po- 
larization of the 21cm radiation is generic for Thomson 
scattering. The cross-correlation will be weaker, for ex- 
ample, if the polarization fluctuations are sourced by 
magnetic fields whose distribution is spatially uncorre- 
cted with the 21cm brightness fluctuations. 

Since Thomson scattering in a nearly uniform medium 
produces E-type polarization and is achromatic, its 
unique fingerprints can be separated from foregrounds 
associated with polarized synchrotron emission. Fara- 
day rotation of the signal as it propagates through the 
IGM will modify the above conclusion, however the fre- 
quency dependence of the effect is precisely known and 
techniques developed to eliminate foregrounds may be 



adapted to reconstruct the original E-type polarization 
signal. It is highly unlikely that any foreground source 
would produce a polarized signal that only contains E- 
type polarization and that correlates with the cosmic 
brightness fluctuations. This may be used as a strin- 
gent test against foreground contamination and instru- 
ment systematics. A similar test has been adopted to 
check if the reconstructed signal from weak gravitational 
lensing surveys is contaminated by instrument al system- 
atics and has proven to be extremely useful (|Refregier I 
l200l . For upcoming experiments which may not have 
the sensitivity to observe the polarized 21cm signal from 
the high redshift IGM, the inferred signal should not cor- 
relate with any observed polarization. This can be used 
to further eliminate the large foregrounds that may con- 
taminate the high redshift signal that is saught after. 
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